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Abstract 

It is widely accepted that topological quantities are useful to describe quantum liquids in low dimensions. The (spin) Hall 
conductances are typical examples. They are expressed by the Chern numbers, which are topological invariants given by 
the Berry connections of the ground states. We present a topological description for the (spin) Hall conductances on a 
discretized Brillouin Zone. At the same time, it is quite efficient in practical numerical calculations for concrete models. We 
demonstrate its validity in a model with quantum phase transitions. Topological changes supplemented with the transition 
is also described in the present lattice formulation. 
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1. Introduction 

Topological quantities are fundamental to describe 
low dimensional quantum liquids where standard sym- 
metry breaking do not have a primal importance [1,2]. 
Typical examples are quantum Hall liquids [3,4,5] and 
anisotropic superconductors with time-reversal sym- 
metry breaking [6,7]. Recently spin Hall conductance 
for semiconductors are also attracting much current in- 
terest [6,7,8,9,10]. The (spin) Hall conductance has a 
characteristic geometrical meaning [3,5]. In some phys- 
ical units, they are given by the first Chern number of 
the Berry connection [4] . 

In practical numerical calculations, we diagonalize 
Hamiltonians on a set of discrete points on the Bril- 
louin zone (BZ). It is thus crucial to develop an effi- 
cient method of revealing the topological property of 
infinite systems with continuum BZ from correspond- 



ing finite systems with discrete BZ. We propose an ef- 
ficient method for the calculation of the Chern num- 
bers on a discretized BZ based on a geometrical for- 
mulation of topological charges in lattice gauge theory 
[11,12,13]. The Chern numbers thus obtained are man- 
ifestly gauge-invariant and integer-valued even for a 
discretized BZ. One can compute the Chern numbers 
using wave functions in any gauge or without specify- 
ing gauge fixing- condittons. Details of the formulation 
and the basic results were published elsewhere [14] . 



2. Topological Description of (Spin) Hall 
conductances 

Chern Numbers as (Spin) Hall conductances: Let us 
consider Chern numbers in the quantum Hall effect as a 
typical example. The spin Hall conductances is treated 
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similarly. We take the BZ by < fe^^ < 27r/gj^j (/^ = 1, 
2 with integers g^). Since the Hamiltonian H(k) is pe- 
riodic in both ki and k2 directions, the BZ is regarded 
as a two-dimensional torus T^. When the Fermi energy 
lies in a gap, the Hall conductance is given by axy = 
— (e^//i) J2n where c„ denotes the Chern number of 
the nth Bloch band, and the sum over n is restricted 
to the bands below the Fermi energy [3]. The Chern 
number assigned to the nth band is defined by c„ = 
^ F with F = dA, where the Berry connection 
A — Ai^dk,^ with Ai^ = {n\d/dk^\n) is defined by a 
normalized wave function of the nth Bloch band |n) 
satisfying H(k)\n{k)) = En{k)\n{k)) and (n|n) = 1 
[3,4,5]. The Chern number can be nonzero only when 
the gauge potential cannot be defined as a global func- 
tion over T^. In this case, one covers by several 
coordinate patches and then, within each patch, one 
can take a local gauge (a phase convention for the 
wave functions) such that the gauge potential is a well 
defined function. In an overlap between two patches, 
gauge potentials defined on each patch are related by 
a U(l) gauge transformation: |n) \n)ij and A — » 
A + uj'^duj {\uj\ = 1). In the continuum theory, one 
needs to fix this gauge freedom to perform any explicit 
calculations. To fix the gauge, one first selects an ar- 
bitrary state which is globally well defined over 
the whole BZ [2] . Then the gauge can be specified by 
|n"*) = P„|0)/Ar'* (if A^'* / 0), where P„ = |n>(n| and 
A'''* = |(<?i|n)| are gauge independent. Generically, this 
gauge is only allowed locally because the overlap N"^ 
may vanishes. It inevitably occurs to have a non van- 
ishing Chern number. Then one needs to use several 
different state \<j)} (gauges) in several regions of patches 
to cover the whole BZ. 

Topological Description on a discretized BZ: We 
now switch to the finite systems with discrete BZ. A 
naive replacement of the difi^erential operator to the 
difference operator breaks the gauge invariance and 
topological characters of the Chern numbers. Here 
we propose au explicitly gauge invariant and topo- 
logical definition of the Chern number on a lattice. 
Let us denote lattice points ke {£ = 1, N1N2) 
on the discrete BZ as ki = (kj-^jkj^), kj^^ = ^-j^, 
(j^ = 0, . . . , — 1). We assume that the state \n{ki)) 
is periodic on the lattice, |n(fc< -I- N^fi)) — \n(ke)), 
where /t is a vector in the direction /i with the magni- 
tude 2n/{q^N^). We first define a U(l) link variable 
firom the wave functions 



Uf,{ke) = N-\ke)(n{ke)\n{ke+fL)}, 

where N^(fc^) = \{n{ki)\n{ki + p,))\. It is well defined 
as long as N^(fc^) 7^ 0. (It is always assumed by an in- 
finitesimal shift of the lattice.) Next wc define a lattice 
field strength by 

Fi2{ke) =\nUi{ke)U2{ke + i)Ui{ke + 2)-''U2{ke)-\ 

(— TT < jFi2{ke) < tt). The field strength is defined 
within the principal branch of the logarithm. Finally, 
we define the Chern number on the lattice as 

We stress here that Cn is manifestly gauge-mvariant 
and strictly an integer for arbitrary lattice spacings. 
To demonstrate it, let us introduce a gauge potential 

Af,{ke) ^lnUf,{ke), 

(— TT < jAij,{ke) < n) which is periodic on the lattice: 
Afi{ke + N^fi) = Aij,{ke). By definition, one finds 

F-i2(ki) =/^iA2(ki) - A2Ai{ke) + 2mni2(ki), 

where A,i is the forward difference operator on the 

lattice, A^f{ke) = f{ki + fi) - f{k(), and ni2(fef) 
is an integer-valued field, which is chosen such that 
(l/i)Fi2{ke) takes a value within the principal branch. 
Now we have 

It shows that the lattice Chern number Cn is an inte- 
ger. To avoid ambiguities we assume that there are no 
exceptional configurations in the system under consid- 
eration, |Fi2(fe£)| ^ TT for all ke [12]. This condition will 
be referred to as admissibility. Under this condition, 
the Chern number is uniquely determined. Since our 
lattice formulation recovers the continuum one in the 
limit — »■ 00, wc expect the lattice field strength F12 
will be small cuougli for a sufficiently large A'^ and the 
lattice Chern number will approach the one in the con- 
tinuum Cn Cn in the this limit. Since both Cn and Cn 
arc integers, we have c„ = c„ for meshes of appropriate 
sizes. As fax as the system is regular, our lattice formu- 
lation could reproduce correct Chern numbers of the 
continuum theory even for a coarsely discretized BZ, 
NlN2^0{\Cn\) [14]. 

Also our method can be extended [14] to the case of 
the non-Abelian Berry connection A = tp^dijj, which 
is a matrix-valued one-form associated with a multi- 
dimensional multiplet ip [2,15]. 
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Uniqueness of the description: ft turns out that the 
description presented so far proves to be unique for dis- 
cretized BZ. Namely, under the admissibility, the space 
of U(l) link variables is divided into disconnected sec- 
tors and the topological number c„ is uniquely assigned 
to each sector. The Chern number c„ is, moreover, a 
unique gauge-invariant topological integer which can 
be assigned to admissible U(l) link variables. The proof 
of this statement has been given by Liischer [11]. 

A key ingredient of his proof is that any U(l) link 
variable can be decomposed uniquely into U^{ke) = 

e'^'^'"^'^'''^+'^'^^''''>A{ke)ujr\ki)A{ke + ft), where 
A^/^"\ke) is a gauge field giving rise to the Chern 
number; Ai A[^"\ke) - Aa^ll^"' (fe^) = 2TviCr,/{NiN2), 
A'^{kt) is a periodic transverse field, A(fc^) is a 
gauge transformation, and Ujr\ki) is a gauge field 
giving nontrivial Wilson lines but giving zero field 
strength. Once A]^"\ki) is found, A'^{ke) can be 
determined uniquely by the relation F\2(ki) = 
Ai (fcf) - A2Aj{ki) + 2iTiCn/{NiN2). Readers who 
are interested in the proof should refer Liischer's paper. 

The admissibility condition is important for the 
present lattice formulation and its breakdown is 
closely related to a singular behavior of the Berry con- 
nection, that is, it is supplemented with a quantum 
phase transition. In the present context of the Berry 
connection, a distribution of the gauge invariant field 
Fi2 is completely governed by the k dependence of the 
Hamiltonian. Each of the topological ordered states 
with a nontrivial Chern number corresponds to non- 
trivial topological sector specified by the admissibility 
and also characterized by the lattice Chern number Cn . 
In the continuum, on the other hand, the topological 
stability of the Chern number is assured by the gap- 
opening condition [2] . The topological quantum phase 
transitions are thus characterized by the gap closing. 
Namely, nontrivial topological sectors of the contin- 
uum, each of which is a topological ordered state, are 
separated by the gaps. 

In passing, we would like to mention the Kubo for- 
mula which is often used in numerical calculations: 

P,^(Jc\ - rHX^ lM"(k)\aiH{k)\m.{k)){m.(k)\a2H{k)\n(k)) 
ri2\l^) — ^'Z>^(^„) [E„(fc)-B„(fe)ia ■ 

This formula is equivalent the field strength F\2{k) 
in the continuum BZ, but not to F\2(ki) in the dis- 
cretized BZ. In the latter case, Fi2{ki) is no longer 
topological, although the gauge-invariance remains. 
Therefore, our topological and gauge-invariant formu- 



lation for discrete BZ should have strong advantages, 
especially of describing regions where topological 
transitions occur. 

3. Quantum Phase Transition and Topological 
Change 

To demonstrate the present scheme, we take the 
Hamiltonian for spinless fermions in an external 
magnetic field: H — tijcle'^'-i cj , where the 

flux per plaquette on the coordinate lattice (j) — 
12n ^ij/i'^^) i'^ P/l- We consider a model with quan- 
tum phase transition, that is, with next nearest 
neighbor (NNN) hopping t' (nearest neighbor hop- 
ping is t). The Hamiltonian in the fc-space is given 
by Hij{k) = ~2t5ijCos{ky - 2n(t)j) — BiSi+i^ — 
B*6i,j+i - B*Si+g-ije~"'''^ -_Bq5,j+g_ie"'''-, where 
Bj =t + 2t' cos [ky - 2'K(i>{j + 1/2)) (i, j = 1, . . . , g) 
with qi — q and g2 = 1 [16]. Bellow, we will present 
some results of applying our method to the middle 
subband of the 0=1/3 (that is, q = 3) system. For 
simplicity, we set A^'i — Nb and A^2 ~ qNs {~ SA^b). 
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Fig. 1. Configuration of integer field ni2(ke) of the NN 
model {t' = 0) in a gauge specified by state \<f)g) over dis- 
cretized Brillouin zones. A'^b = 3(a), 4(6), 5(c) and ll(ci). 
Black (white) circles denote ni2 = —1 (1). 

We show the integer field ni2{ke) in Figs. l(a)-(d), 
where we use the global gauge specified by the states 
|<j!>9) = e'''(*-+'=«'(l, 1, 0)"^. The black and white circles 
denote ni2 = — 1 and 1, respectively, whereas a blank 
implies ni2 = 0. It is clear that any of them gives the 
correct Chern number c„ = —2. The field ni2(ki) is 
gauge-dependent, but their sum is gauge-invariant. It 
also shows a convergence of the gauge dependent field 
ni2 in the limit Ab — » oo. Although we have fixed 
the gauge to calculate the field ni2 , calculations of the 
Chern numbers are performed in any gauge. We do not 
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need specific gauge-fixing to make the gauge connec- 
tion smooth. An arbitrary gauge (e.g., a phase choice 
of eigenvectors given by a numerical hbrary) can be 
also adopted to compute the Chern number. 

Now let us change a ratio t' /t in the NNN model. In 
Figs. 2, we show the energy spectra of the NNN mod- 
els. They are modified Hofstadter's butterfly diagrams. 
The integers in the figures are the Hall conductances 
when the fermi energy lies in the gap, that is, it is a 
sum of the Chern numbers below the energy gap. Due 
to the sum rule of the Chern number [2], the Chern 
number of the specific energy band is given by a differ- 
ence of two integers, above and below the band. As is 
known [16], the NNN model shows a topological quan- 
tum phase transition accompanying a discrete change 
of the Chern numbers. For example, an increase of t' /t 
causes a rearrangement of the gap connectivity near 
t' /t — 0.26 at = 1/3. One can see that the energy gap 
just above the second band at = 1/3 is rearranged by 
the small change of t' /t. It is associated with a change 
of the gap labeling from —1 to +2. Since the large en- 
ergy gap just below the middle band is stable and the 
gap is labeled as -|-1, the Chern number of the middle 
band changes from —2 to -|-1. This topological change 
of the Chern number is confirmed by a calculation of 
the gauge dependent field ni2 in Figs. 3(b) and (d). The 
field strength F12 just before and the after the quan- 
tum transition in also shown in (a) and (c). A small 
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Fig. 3. (a):latticc field strength F\2 and (b): configuration 
of integer field n\i{ki) of the NNN model with t' = 0.267t, 
ATg = 15. (c) and (d) the same for parameters t' = 0.268t. 

model ) occurs at t' = tc (0.267 < t'Jt < 0.268). This 
quantum phase transition is associated with a topolog- 
ical change of the integer value field ni2 as shown. 

This work was supported in part by Grant-in- Aid for 
Scientific Research from JSPS. 



Fig. 2. Energy spectra of the NNN models as a function of 
the flux per plaqucttc (p. (a) t' = O.lSt and (b) t' = 0.32t. 
Rearrangement of the gap connectivity is observed at the 
second gap near <p = 1/3. 

change of the parameter t' /t does not change the spec- 
trum so much. On the other hand, the field strength 
F12 is singular at the critical point. It is consistent with 
the breakdown of the admissibility as expected. That 
is, the quantum phase transition ( of the present lattice 
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